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We prove that static, spherically symmetric, asymptotically flat, 
regular solutions of the Einstein- Yang- Mills equations are unsta- 
ble for arbitrary gauge groups, at least for the "generic" case. 
This conclusion is derived without explicit knowledge of the pos- 
sible equilibrium solutions. 
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1 Introduction 

In several recent papers ||l|, 0, ^ |^ we have studied important aspects of the 
Einstein- Yang-Mills (EYM) system for arbitrary gauge groups. In particu- 
lar, we investigated the classification and properties of spherically symmetric 
EYM solitons (magnetic structure, Chern-Simons numbers) and a general- 
ization of the Birkhoff theorem for the non-Abelian case. We also worked 
out the generalization of the first law of black hole physics (Bardeen-Carter- 
Hawking formula), allowing for additional Higgs and dilaton fields |^ |^. For 
other studies of these and related topics we refer to [^ ||, 3 . 

In the present paper we prove that static, spherically symmetric, asymp- 
totically flat, regular solutions of the EYM equations for any gauge group 
are unstable, at least in the "generic" case (defined below). In a recent letter 
[ p!0|] we have already sketched how we arrived at this result. Here, we present 
full details of the proof and discuss also some further mathematical issues 
involved. 

On physical grounds, this instability was expected, because we had shown 
earlier that the Bartnik-McKinnon solutions |TI| for the gauge group SU(2), 



as well as the related black hole solutions [|12], |T^, |T4[ are unstable |T^, |T6|, 
P!7| , |TB| . A mathematical proof of this expectation presents, however, quite a 
challenge, since one can not rely on any knowledge of the possible solutions 
(apart from regularity and boundary conditions). 

Our strategy is based on the study of the pulsation equations describing 
linear radial perturbations of the equilibrium solutions and involves the fol- 
lowing main steps: First, we show that the frequency spectrum of a class of 
radial perturbations is determined by a coupled system of radial "Schrodinger 
equations". Eigenstates with negative eigenvalues correspond to exponen- 
tially growing modes. Using the variational principle for the ground state it 
is then proven that there always exist unstable modes (at least for "generic" 
solitons) . 

There is, unfortunately, no direct way to apply our method to black 
holes, because of problems related to the boundary conditions at the horizon. 
We have, however, recently used a similar procedure to establish also the 
instability of gravitating regular sphaleron solutions of the SU(2) Einstein- 



Yang-Mills-Higgs system with a SU(2) Higgs doublet |19], which have been 
constructed numerically in 



The paper is organized as follows: In Sec. 2 we recall some basic facts 
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and equations of our previous work 0, Q which will be needed in the present 
analysis. In Sec. 3 we then derive the linearized perturbation equations 
for solitons (and black holes) and bring them into a convenient partially 
decoupled form. The resulting eigenvalue problem is discussed in Sec. 4, and 
in Sec. 5 we show the existence of unstable modes. 

2 Spherically symmetric EYM fields 

We begin with a convenient description of gauge fields with spherical symme- 
try (for derivations see 0). Let us fix a maximal torus T of the gauge group 
G with the corresponding integral lattice / ( = kernel of the exponential map 
restricted to the Lie algebra LT of the torus T). In addition, we choose a 
basis S of the root system R of real roots. The corresponding fundamental 
Weyl chamber 

K{S) = {H eLT\ a{H) > for all a G 5 } (1) 

plays an important role in what follows. We have shown in [Q that to a given 
spherically symmetric gauge field configuration, there belongs in a natural 
way a canonical element Hx G lr\K{S) which characterizes the corresponding 
principal bundle P{M, G) over the spacetime manifold M admitting an SU(2) 
action. If the solution is regular at the origin, H\ lies in a small finite subset 
of / n K{S), which we have described in 0. In much of our discussion we 
exclude (for technical reasons) the possibility that Hx lies on the boundary 
of the fundamental Weyl chamber. The term generic always refers to fields 
for which the classifying element Hx is contained in the open Weyl chamber 
K{S). 

The SU(2) action on P{M, G) by bundle automorphisms induces an ac- 
tion on the base manifold M. An SU(2)-invariant connection in P{M, G) 
defines a connection in each subbundle over a single orbit of this induced ac- 
tion, which by Wangs theorem is described by a linear map A: LSU(2) -^ LG, 
depending smoothly on the orbit and satisfying 

Ai = [A2, A3] , A2 = [A3, Ai] , A3 = ~Hx/47c . (2) 

Here, A^ := A(r,fc) and 2iTk are the Pauli matrices. These equations imply 
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that A+ := Ai + 2A2 lies in the following direct sum of root spaces La of LGc'- 

^+ £ ^ La , 

oGS (3) 

E:={aeR+\ a{Hx) = 2}. 

R+ denotes the set of positive roots in R (relative to the basis S). In the 
generic case S turns out to be a basis of a root system contained in R. (This 
is proven in Appendix A of Ref. 0].) 

The LG-valued functions on the orbit space determine part of the con- 
nection on P{M, G). Before giving a parametrization of the YM fields in a 
convenient gauge, we fix our conventions in parametrizing the Lorentz metric 
g on M which is, of course, assumed to be invariant under the induced SU(2) 
action. We use standard Schwarzschild-like coordinates and set 

g = -NS'^df + N-^dr^ + r2(t/^2 ^ ^^^^2^ ^^2^ ^ ^4^ 

where the metric functions N and S depend only on r and t. We use also 
the usual mass fraction m{r, t), defined by A^ =: 1 — 2m/ r. 

In Ref. it is shown that there exists always a (local) gauge such that 
the gauge potential A takes the form 

A = A + A, (5) 

with 

i = A2rf^+ (Agcos^- Aisin^)rfv? (6) 

and 

A = NSAdt + Bdr, (7) 

where A and B commute with Hx (i.e. with A3). If H\ is generic its central- 
izer is the infinitesimal torus LT. Hence, A and B are LT-valued and A is 
thus Abelian. 

For the example of the gauge group SU(2), H\ is an integer multiple of 
47r T^: Hx = 4:7ik T3 with k E Z, and the only solutions of (^ are Ai = A2 = 0, 
A3 = fcrs, or 

Ai = WT1+WT2 , A2 = TW Ti±WT2 , A3 = ±T3 . (8) 

We introduce some further notation which is used below. A suitably 
normalized Ad(G')-invariant scalar product on LG will be denoted by (•,■). 
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We use the same symbol for the hermetian extension to LGc (antihnear in 
the first argument), and | • | means the corresponding norm. Note that the 
original Ad((j')-invariance extends on LGc fo 



■X,[Z,Y]) + {[c{Z),X],Y) = 



(9) 



where c is the conjugation in LGc- 

If we insert the parametrizations (§)-(0) into the EYM equations, we 
obtain a system of partial differential equations for the metric functions A^, 
S and the YM amplitudes A±, A, B. For our instability proof it will suffice 
to write them in the temporal gauge .4 = 0. Specializing the results of pf 
(and using slightly different notation) they read as follows: 

The Einstein equations give two constraint equations for the r derivative 
(denoted by a dash) and the t derivative (denoted by a dot) of m 



m 



{NG + r^pe}, 



K 



m 



-NH 



K :-- 



8nG), and the (rr )-equation reduces to 



S 



-G 



where 



G 
H 

Pe 



Re(A+,AV+[S,A+]) 

— (l^l' + l^lH 



[S,A_ 



with 



^ 



E 



T 



S" ' ' 2 
The YM equations decompose into 



[A+,A- 



A, 



(10) 

(11) 

(12) 
(13) 

(14) 

(15) 



'■r-B) + 



NS\S 



A+ , A'_ + [S, A_] + A_,AV + [S,A 



, (16) 
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1 / 1 : 

S\NS ' 



Uns{a'^ + [b,a^]})' 



S 

Ni^[B,A'^] + [B,[B,A 

2 1 



2r-B)^2^-[B,B,^- 



+ -[^,A+] = 0, (17) 



[A+,A_] + [A_,Ah 



0. 



:i8) 



The last equation is the Gauss constraint. For the generic case the term 
proportional to [i3, i3] in ([T8|) vanishes. 

For static solutions all time derivatives disappear and the basic equations 
simplify considerably. (For the Bartnik-McKinnon solutions A is of the form 
(D with w = 0, A3 = Tg and i = in (§).) 

3 Perturbation equations 

In this section we study time-dependent perturbations of a given static, 
asymptotically fiat solution of the coupled EYM equations (|T0D , ([TT|) , (|16])- 
(|TBp. Regular solutions are of purely magnetic type (^ = in (^) with 
vanishing YM charge. Unfortunately, this is not yet rigorously proven under 
satisfactory weak fall-off conditions, but there is strong evidence for this (see 
[^, |21| for partial results.) The perturbation equations which we shall derive 
hold also for black holes, if these are assumed to be of purely magnetic type. 

From now on the symbols A±, N, S, etc. refer to the equilibrium solution 
and time-dependent perturbations are denoted by 6A±,6B, etc.. All basic 
equations are linearized about the equilibrium solution. 

First, we linearize the right hand sides of the Einstein equations ([T0|) and 
(|TT]). Since B and A± vanish for the equilibrium solution, the first order 
variation of the source G is 



5G= Re(A'^, <5A'+)- Re(A;, [A+,5B]) 



(19) 



Here, the last term vanishes, because the properties of the scalar product 
mentioned earlier (notably (01)) imply that 



-2Re(A;, [A+,5S]) = ( [A+,A'_] + [A_,A'J,^-S) 



(20) 
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and the YM equation (|T6|) for the equihbrium solution shows that 

[A+,A'_] + [A_,AV] = 0. (21) 

Thus 

SG= Re ( A'^ , 6A'^ ) . (22) 

The only first order variation for pg conies from S\J-'\ = 2{j-' , SJ-"), with 
(see eq. (|15D) 

6:F='-[A+,6A.]-'-[A.,6A^], (23) 

giving 

6pe = ^^Re{t[T,A+],6A+). (24) 

Now we can work out the variation of the first Einstein equation in ([T0|) . 
With (|22|), ( p^ and (0) for the equilibrium solution, we find 

6m' = -j6m+'^^NRe{A'^, 6A'^) + Re{ -^ [J^, A+] , (5A+ ) | . (25) 

For the commutator in the last term we use the unperturbed YM equation 
fUl, i.e. 



or 



^ [:F, A+] = NjA'_, + N'A'^ + NAl , (26) 
whence 

W = -|(5m + |{|iVRe(A^,M+)} + |^ArRe(A;,<5A+)}' (27) 

{6mSy = (-NSRe{A'^,6A+)X . (28) 
Therefore, Sm must be of the form 

6m = ^NRe{A'^,6A+) + ^, (29) 

where f{t) is a function of t alone. This function is determined by considering 
the variation of the second Einstein equation in (0), which reads 

6m='^NRe{A'^,6A+). (30) 
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Sm='^NRe{A'^,SA+)+g{r), (31) 



with a function g{r) of r alone. By comparing (29) and (31) we arrive at the 
remarkably simple result 

Sm=-NRe{A'^,6A+) , (32) 

which generalizes an observation already made in |jT5[ . 



The variation of the Einstein equation (|Tl|) is immediately obtained with 

m 

<^(§)=^iVRe(A;,MV). (33) 

Before also linerizing the YM equations we introduce a decomposition of 
A+ and (5A+ into "real" and "imaginary" parts. For this we introduce (for 
a G S) a basis e^ of the root spaces La in the direct sum (|^) with respect to 
which we expand the unperturbed A+ as well as its perturbation (5A+, 



A+ = ^ w°ea, SA+ = Y, 5w'^ e„ . 


(34) 


Then we have 

5A± = (5X± ± i5Y^ 


(35) 


with 

5X+ = "^ Re {5w'^) Ca , 5Y+ = Y, Im {5w'') e„ 


(36) 



and the corresponding expansion for 5X_ and 5Y- with Cq, replaced by c{ea) G 
L_Q,, because 5A_ = c((5A+) and thus 

6X_ = c{6X+) , 6Y_ = c{6Y+) . (37) 

We shall call 6X±, 6Y± the real and imaginary parts of the perturbations 
6A±. It was shown in 0] that the unperturbed A+ can be chosen to have 
only a real part. 
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This decomposition will lead to a significant decoupling of the perturba- 
tion equations. Note in particular, that the variations 6m and Spo in (|^) 
and (^) depend only on the real part 6X+: 

Sm = |iV(A;,5X+), (38) 

Spe = -Az[:F,A+],6X+). (39) 



We consider now the first variation of the YM equation (0). Its decom- 
position into real and imaginary parts yields 

^ 5X+ = -N6XI- ^^6X'^- ^[A+,6:F] + ^[:F,6A+] 






6NA'[ - 6(^^^)a'^ (40) 



and 

~ NS 



L^ 5f+ = -iv| 5F; + t [A+, 6B]' + t [AV, 6B] 

-^-^[6Yi + ^[A^,5B]^ + ^^[:F,5Y^]. (41) 

The third term on the right hand side of (|40| ) is indeed real and can be 
written, using (^2]), as 

^ [A+, 5J^] = 4 ad(A+) ad(A_) SX+ . (42) 



Equation (|40|) can be simplified further. From (|38D and the equilibrium 



equation (^) we deduce 



- 5N A'i = -5mA'' 
r 



fi;A^Re(A'+, 5X+) A"_^ 
/.Re(AV,<5X+)|-^A+ + -^[^,A+]} 
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and the Einstein equations (p!OD, ( p^ give 



5 ^ = — :^6m + Kr 5pe • (43) 



If we use here (|38D and (|39D we see that the last two terms in (^OD can be 
expressed as follows: 



y-(-^^)7{^+^}<-^-^^^ 



_ SNAl - 5^ = ^^<i -(P.A+) 'ii ^^^ + ^\{ p.K , 5X^ ) 



+ (P.A+) ^( [^, A+] , 5X+ ) + [P, A+] ^( j9,A+ , <5X+ ) I , (44) 
where we have introduced the differential operator 



Inserting these expressions into (^0]) gives the following pulsation equation 
for the real amplitude 6X^ of the YM field 

6X+ + Uxx SX+ = 0, (46) 

where the operator Uxx is given by 

Uxx = p.' + ^;^ad(^^) - ^^ad(X5A+)ad(X5A_) 

+ (p.A+)^([^,A+],-) + [^,AJ^(p.A+,-). (47) 

It is remarkable that the perturbations 6Y± and 63 do not appear in (^) 
and that the back reaction of gravitation on 6X^ can be described by an 
effective potential (last three terms in (|i7|)). 



Equation (^T]) can easily be brought into the form 

6Y+ + Uyy SY+ + UybSB = 0, (48) 
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where 



Uyy 

Uyb 



2 NS^ ,,.^, 
p^ad{NSA+) + ad{NSp^A- 



(49) 
(50) 



We have thus achieved a partial decouphng, because neither 5X+, nor the 
metric perturbations appear in (|48D. 

We proceed with the hnearization of the YM equation (p^ . The variation 
of the last two terms is 

- "a+, [A-,6B]] + [A+,(5A'_] - [A'_,5A+] + conjugate, (51) 

which leads (with 6A± = 6X± ± iSY±) to 

A+ , [A_, 6B]]+t [A+, 6YL] + % [A'_, 5r+ 



+ I [A+, bX'_\ - [A'„, (5X+] I + conjugate . 



Here, the terms in the first curly bracket are in LT ^ while those in the second 
are in iLT . The latter are compensated by their conjugates and we find 



Nr'^m^ + Ubb SB+ + Uby 6Y+ = 



with 



Use = - ad(iV5A+) ad(A^5A„) , 

Uby = -ad{NSA_)p, + ad{NSp,A_) 

At this point we collect the results obtained so far as follows: Let 



(52) 



(53) 
(54) 



5* 




T 




(55) 



then (^) and (|S2|) can be written as a 2 x 2 matrix equation 

T6^ + U6^ = , 



(56) 
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with 

U=\ ^^^ ^^^ I . (57) 

Uby Ubb 



The operators in this matrix are given in eqs. (|49|), (|50|), ([53|) and (^41). 

The perturbation equations (^) and ( |5BD do not include the Gauss con- 
straint (|T8|), whose hnearization is easily found to be 

^|p.(^5i3)-[A+,5r_]}=0. (58) 

The role of this constraint will be discussed later. 

In concluding this section we emphasize once more, that the perturbation 
equations hold also for black holes, if these are assumed to be of purely 
magnetic type. In our further discussion we will, however, consider only 
perturbations of uncharged regular solutions. 



4 The eigenvalue problem 

It is natural to introduce the following scalar product for LGc-valued func- 
tions on IZj^: 

(01^)= / (0,^) — , (59) 



/o '" "' NS 

because the operators Uxx^ U and T are symmetric with respect to this 
scalar product on a dense domain of L^-functions. This can easily be seen 
by using 

(0b*^) = (p*0IV^) (60) 

for smooth functions which vanish at the origin, and 

(0|ad(Z)V^) = -(ad(c(Z))0|^) (61) 

for all LGc-valued functions 0, ?/^, Z in L^ (see (^). 
We specialize now to harmonic perturbations 

5X+(r, t) = ^(r)e~^'^* , (5*(r, t) = 6^{r)e-'^' , (62) 

whose frequencies satisfy the eigenvalue equations 

Uxx^ = oj\ , (63) 
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and 

US^ = uj^T6<^. (64) 

It should be remarked at this point that P7| ) and (^) reduce for the 



Bartnik-McKinnon solution [|rT| to the eigenvalue problem derived in Ref. 



ITH] , where it was shown that this has exactly one unstable mode. (A similar 



instability for the "colored" black hole was found in ||16|| .) 

Let us now turn to the role of the linearized Gauss constraint (^) in 
conjunction with the eigenvalue problem (0). We show first that a variation 
(5$ is orthogonal with respect to the scalar product defined by T, 

(■|-)t=(-|T|.), (65) 

to all gauge variations 

^'^gauge = ' , ^ ' '^^^^ 

if and only if the curly bracket in (|58| ) vanishes. Note that these gauge 
variations arise if (||) is subjected to the gauge transformation g = exp(ex), 
because (^ and (|^ show that this induces the infinitesimal transformation 

A+^A+-[x,A+], B^B + x'. (67) 

To prove the statement just made, we compute 

/■CO , ^ fjf 

(5$|r|<5$gauge) = I {{5Y^,-t[A^,x]) + {SB,Nr\')} — , 

Here we have used (|^) and made a partial integration, dropping a boundary 
term. This is allowed if x is regular at the origin and vanishes sufficiently 
fast at infinity. Since x is otherwise arbitrary and i [A_, 6Y_^] = —i [A_|_, 5Y_], 
eq. (l68| ) implies our claim. 

In the next section we will show that the eigenvalue equation (^4[) has 
at least one mode 5<^ with cj^ < 0. Such a mode is orthogonal with respect 
to the scalar product (|65| ) to any zero mode of ( |64D and thus in particular 
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to (5$gauge in (|66D. (This follows since different eigenmodes of (p^ ) are 
orthogonal with respect to the scalar product (|65|) , because U and T are 
symmetric with respect to (§^).) Hence, we can conclude that the Gauss 
constraint is automatically satisfied. 
From (|6^) we obtain 

and for the frequence Uo of the fundamental mode we have the minimum 
principle 

We do not discuss here the precise mathematical nature of the eigenvalue 
problem (domains of definition, essential selfadjointness, etc.), because the 
functional analytic aspects are very similar to other well-studied eigenvalue 
problems. 

5 Instability of generic EYM solitons 

We are now ready to establish the main point of this paper: 

For a given regular solution with A+ = X^ogs W^^a we shall construct a 
one-parameter family of field configurations A(^)_|_, i3(^) such that the vari- 
ational expressions ([70|) for 6A± = {dA(^^)±/dx)x=o, ^^ = {dB(^x)/dx)x=o is 
negative . This family is chosen of the following form: 

A(^)+ = Ad(exp(-x^)){A+cos(x)+^A+(oo)sin(x)}, (71) 

%) = XZ', (72) 
where Z is an LT-valued function of r with the properties 

lim [Z,AJ=iA+{oo) , suppZ' C [1 - e , 1 + e] (73) 

r— >0,oo 

for an e > 0. The existence of such a function can be seen as follows: Let 
{ha}aeT: be the dual basis of 27rS and put 

y-X^yau 7a _ / w°(oo)/w"(0) forr<l-e, . . 

^-^^^ ^"' ^ "t 1 forr>H-6. ^^^^ 
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It is easy to verify that both conditions in (|73D are satisfied. (In Appendix 
A of ref. [Q we have shown that w"(0) 7^ for all a G S.) 

We note a few properties of the family (^), (^). The equilibrium solution 
is clearly obtained for x = 0. Applying a gauge transformation with g = 
exp(— x^) we obtain 

^{x)+ -^ ^+ cos(x) + iA+{oo) sin(x) , %) -^ . (75) 

The first variations of (|7TD and (^) are 

M+ = -[Z, A+] + ?A+(oo) , 6B = Z' , (76) 

and these satisfy by construction the desired boundary conditions 

lim (5A+ = 0, lim 6B = . (77) 

r^OjOO r— >0,oo 

{6B vanishes even outside [ 1 — e , 1 + e ] .) Finally, (5A+ has only an imaginary 
component 

6Y+ = -i5A+ = i [Z, A+] + A+(oo) (78) 

and thus by ([70|) 

with (5$ = ((5K|., 5S) given by ( ffSD and the second eq. in (176|). 

This judicious choice of trial functions fulfills our goal: The denominator 
in (^) is finite and the numerator turns out to be strictly negative! 

The first of these two points is simple. Since 6B in ([76|) has compact 
support, we have to check only whether 

{6Y^,SY^)— < 00. (80) 

By construction, i6Y+ = A+(r) — A+(oo) for r > 1 + e. Since A^ and S both 
approach 1 at infinity, the integral is finite if A_|_(r) — A+(oo) is assumed to 
converge to zero faster than r^^^"^. 

The calculation of the numerator in ( |79D is somewhat tedious. Consider- 
able simplifications occur by separating a gauge mode in 5$: 

/ A+(oo) \ / i\Z,A+i , 

6^ = d^$gauge + , ^-^gauge = | | • (81) 
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Clearly U 6^ gauge = 0, and thus (|9|) and d^) give 

f/yyA+(oO 



I TL...\ (,^^\ \ AT02 I 



U5<^ 



From (^) and (82) we obtain 



{6<^\U\ 6'^ ) 



_NS^l z[^,A+(oo)] 



+ / (A+(oo),2[^,A+(oo)] )- rfr 



+ f {Z',-i[A'_,A+{oo)])NS dr 



(82) 



(83) 



Let us denote the integrands of the three terms by Ji, J2 and J3. We find 
immediately 



J^ = -(^,[A+(oo),[Z,A_]]) 
J2 = 2^(^,A3). 



(84) 
(85) 



In the second equation we have used (|T5|) and the vanishing of the YM 
charge, implying that lim^^oo A(r) is a homomorphism from LSU(2) to LG, 
whence 

i[A+(cx)),A_(oo)] = 2A3. (86) 

In a next step we show that the first and the last term in ( |5BD compen- 
sate each other. For this we rewrite the third term, performing a partial 
integration and making use of the equilibrium equation (^61) , as follows 



J-^dr 



i / {NSA'^, [Z,A+(oo)]') dr 



t{NSA'^, [^,A+(oo)]) 



00 




S 



(^,[A_, [Z,A+(oo)]])-dr . (87) 
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The boundary term vanishes and the double commutator in the last term is 
equal to [A+(cx)) , [-Z', A_] ], as is seen by using the Jacobi identity and the 
fact that [A+(oo), A_] is in iLT. Comparing this result with ( P^ ) shows that 
there remains indeed only the second term in (|83|). Thus from (|85|) we obtain 
the intermediate result 

S 



'6<^\U\6<!>) = 2 f {:F,A3) — dr. (88) 



Finally, we show that the last expression can be transformed into a form 
with a definite sign: 

2 J CF , A,) ^^dr = - J {NS\A'^f + 2^^\:Ff} dr . (89) 

In order to see this we perform a partial integration of the first term on the 
right and use again the YM equation (^): 

l{A'^,NSA'^)dr={A+,NSA'^)\^ - J {A+ , t[T,A^] ) ^ dr . (90) 

The boundary term vanishes and the integral combines with the last term in 
(|89D to the left hand side, because we have (see (|15|)) 

2\:Ff = {:F, I [A+,A_] - A3) = (A+, z [^, A+] ) - (^, A3) . (91) 

All together, we have established the crucial result 

(5$|t/|(5$) = - I [nS\A'^\^ + 2^\:F\^] dr . (92) 

This expression is finite and strictly negative. Hence we have shown that ujq 
is indeed negative, and thus there exist unstable modes of (|6^) . These fulfill, 
we recall, automatically the linearized Gauss constraint (pS]). 

One can show that the expression ( p2|) is also equal to the second variation 
of the Schwarzschild mass for the one-parameter family (p9|), (|72|). (This 



is the way we arrived originally at the variation (|72D ). For a systematic 
discussion of the relation between variational principles for the spectra of 



radial pulsations and second variations of the total mass, we refer to |22 . 
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In summary, we have proven the instabihty of all generic, regular equi- 
librium solutions. More precisely, we have established the following 

Theorem. A static, sperically symmetric, asym,ptotically flat, regular solu- 
tion of the EYM eqs. ([Tdj), ^Tlp, (J7^j-(j7^j is unstable if the following three 
conditions are satisfied: 

(i) The solution is generic (i.e. the classifying element H\ = — 47rA3 lies 
in the open fundamental Weyl chamber). 

(a) The (magnetic) YM charge vanishes (i.e. linir^oo A(r) is a homomor- 
phism from LSU(2) to LG). 

(Hi) Asymptotically A+(r) — A_|_(oo) ~ r~" with a > 1/2. 

We emphasize again the strong evidence, that the assumptions of the 
theorem together with condition (i) already imply condition (ii). Also, we 
would like to stress that we were able to draw this conclusion assuming 
only weak asymptotic conditions for the solitons. In particular, the fall-off 
conditon (iii) of the theorem is mild and is certainly fulfilled for the Bartnik- 
McKinnon solutions, as was shown rigorously in [^. The same is true for the 
solutions which have been found numerically by H.P. Kiinzle for the group 
SU(3) [pi|]. (For both types the exponent a is equal to one.) 
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